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Abstract 

We show that in the presence of sufficiently strong v T annihilations into ma- 
jorons, the primordial nucleosynthesis constraints can not rule out v T masses 
[ in the MeV range. 
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1 Introduction 



Despite great experimental efforts the tau-neutrino still remains as the only one that 
can have a mass in the MeV range. The present experimental limit on its mass is 
1 

m Vr < 23 MeV (1) 

Further progress will have to wait for the tau-charm or B-factories. On the other 
hand, many particle physics models of massive neutrinos lead to a tau neutrino with 
mass in the MeV range. Moreover such a neutrino may have interesting cosmological 
implications. It is therefore important to examine in more detail the cosmological 
constraints. 

The first comes from the critical density argument. From the contribution of 
stable v r to the present relic density one gets m Vr < 92Qh 2 eV, where h is related 
to the Hubble constant through h = /Zo/(100 kms~ l Mpc~ l ). This means that a 
massive v T with mass in the MeV range must be unstable with lifetimes smaller than 
the age of the Universe. However, it has been shown that in many particle physics 
models where neutrinos acquire their mass by the spontaneous breaking of a global 
lepton number symmetry^] this limit can be avoided due to the existence of fast v T 
decays. 

The second constraint comes from primordial nucleosynthesis considerations. In 
the standard model, these rule out v T masses in the range |J 

0.5 MeV < m Ur < 35 MeV (2) 

It is possible to weaken the constraints on u T of Eq. @, by adding new interactions 
beyond the standard ones. One may either consider that the z/ T 's are unstable 
during nucleosynthesis || or that they possess new channels of annihilation beyond 
the standard ones. It is this last possibility that we investigate in the case where 
the can annihilate to majorons (J). 

2 Evolution of v T number density 
2.1 Before weak decoupling 

We will assume that the massive Majorana u T y s are stable during the nucleosynthe- 
sis epoch. They interact with leptons via the standard weak interactions u T T\ «-> 
vqVo, e + e~. Moreover the z/ r 's annihilate to majorons v T v T — > J J via the diagrams 
shown in Fig. (H). The t-channel diagram is present in all Majoron models, while the 
strength of the s-channel scalar exchange diagram is somewhat model-dependent. 
For the range of g values relevant for our purposes (g > 10 -5 ) it is quite reasonable 
to neglect the s-channel process, especially in Majoron models where the breakdown 
of the global lepton number symmetry happens at the electroweak scale or higher, 
like the simplest seesaw model JJ. We have shown [J7| that for models with low 
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Figure 1: Feynman diagrams for annihilations of tau neutrinos into Majorons. 

scale breaking of lepton number, such as given in ref. || [| the bounds obtained 
by neglecting the s-channel process could be relaxed by at most a factor ~ 2 or 
so. Therefore in what follows we will consider that the z/ r 's annihilate to majorons 
v T v T — > J J via the diagonal coupling 



C VtVtJ = i\g J Vrlh^r 



(3) 



The evolution of the number densities of v T (n T ) and majorons (nj) is given by 
the solution of a set of Boltzmann differential equations 



h T + 3Hn T = - £ ((Tiv) (n 2 T - «?) 2 ) - (a jV ) (n 2 T - «") 2 - 
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where {&iv) are the thermally averaged cross sections times the v T relative velocity 
v. The cross sections are for annihilation to majorons 
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where b 2 L + b 2 R = 1/2 for i = u and b 2 L + b 2 R = 2 [(-1/2 + sin 2 6 W ) 2 + (sin 2 6 
0.25 for i = e. In Eq. ([|) we assumed Boltzmann statistics and rij = for i = e, u . 

Now let us describe our calculations. First we normalized the number densities to 
the number density of massless neutrinos, n ~ 0.181T 3 , introducing the quantities 



2 



r a = n a /no, where a = u T , J, and the corresponding equilibrium functions r^. Then 
from Eq. (||) we get 

^ = (^-3Hr a ) i-|r Q (8) 
dT \n J T T K J 

On the other hand, the time derivative of the temperature is obtained from Einstein's 
equation 

p = -3H(p + P) (9) 
where p is the total energy density and P is the pressure. 



2.2 Past weak decoupling 

Once the z/ T 's decouple from the standard weak interactions, they remain in contact 
only with the majorons. Then one has two different plasmas, one formed by z/ T 's and 
J's and the other by the rest of particles, each with its own temperature, denoted 
by T and T 7 , respectively. We assume that the photon temperature evolves in the 
usual way, y = Hy. The evolution equations of the v T and J number densities are 
now simplified versions of Eq. (||), because we S Ut = —Sj, 

h T + 3Hn T = —Sj 

hj + 3Hnj = Sj (10) 

To solve these equations one must find a relation between T and T 7 . This is obtained 
using Einstein's equation for the v T + J plasma. 

In order to determine the final v T frozen density which will be relevant during 
nucleosynthesis we have to solve numerically the appropriate set of coupled differ- 
ential equations for each pair of values (m UT , g) . Before weak decoupling these are 
Eqs. (|) e (|) while after decoupling one should use Eq. ([[(]). The initial conditions 
are, for sufficiently high temperatures, r Ur = r^, rj = r e j and T Vt = Tj = T uo . 



3 Nucleosynthesis constraints 

The value of r UT (m UT ,g) is used to estimate the variation of the total energy density 
Ptot = Pr + Pu T - In Pr ah relativistic species are taken in account, including the ma- 
jorons and two massless neutrinos, whereas p VT is the energy density of the massive 
z/ r 's. 

In order to compare with the standard model situation, we can now express 
the effect of the v T mass and of the annihilation to majorons in terms of an effec- 
tive number of massless neutrino species, N eq , which we calculate for each value of 
r v T ( m u T ,g)- To do this, we first numerically calculate the evolution of the neutron 
fraction, r n , by varying the value of iV e g[[lO]. Then we incorporate ptot in this numer- 



ical code and perform the integration for each pair of (m UT ,g) values. Comparing 
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Figure 2: Effective number of massless neutrinos equivalent to the contribution of 
heavy zz/s with different values of g in units of 10~ 5 . For comparison, the dashed 
line corresponds to the standard model case when g = 0. 

the r n obtained in each case at T 7 ~ 0.065 MeV (the moment when practically all 
neutrons are wound up in 4 He), we can relate (m„ r , g) and N eq . 

The results are shown in Fig. (|2]). One can see that for fixed N eq , a wide range 
of tau neutrino masses is allowed for large enough coupling constants. In fact all 
masses below 23 MeV are allowed provided the coupling constant g exceeds a value a 
few times 10~ 4 . For comparison, the dashed line corresponds to the case when g = 
and no majorons are present. These results can also be expressed in the m Vr — g 
plane, as shown in Fig. @. 




4 Majoron models 



There has been a variety of majoron models proposed in the liter at ure|TT|. They 
are attractive extensions of the standard model where neutrinos acquire mass by 
virtue of the spontaneous violation of a global lepton number symmetry. Apart 
from a phenomenological interest of their own, majoron models offer the possibility 
of loosening the cosmological limits on the neutrino masses, either because neutrinos 
decay or annihilate to majorons. 

We have seen above that the restrictions imposed by primordial nucleosynthesis 
upon a heavy tau neutrino disappear for sufficiently large values of the v T v T majoron 
coupling g. Different models imply different expectations for the coupling g and for 
the relation between g and the v T mass. Just to give a concrete example let us 
consider the supersymmetric models with spontaneous violation of R parity In 
these models we get the relation between g and m Vr depicted in Fig. (HI). The 
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Figure 3: The values of g(m UT ) above each line would be allowed by nucleosynthesis 
if one adopts the N™ ax = 3, 3.4, 3.8, 4.2 (from top to bottom). 

different curves correspond to different values of the parameter vr of the model. We 
see that we can obtain for g values of the order of a few times 10~ 4 , as required in 
our nucleosynthesis analysis. 

5 Conclusions 

We have determined the restrictions imposed by primordial nucleosynthesis upon 
a heavy tau neutrino in the presence of sufficiently strong v T annihilations into 
majorons. We have shown that for values of the v T v T majoron coupling in excess 
of 10 -4 all the masses up to the experimental limit of 23 MeV can be allowed. We 
have shown that in a supersymmetric model with spontaneous violation of R parity 
such values can be obtained. 
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